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Abstract 

We show within a geometrical model developed in earlier papers that mul- 
tiplicity distributions are cut off at large multiplicities. The position and 
motion of the cut-off point is related to geometrical- and KNO scaling and 
their violation, in particular by the rise of the ratio a e i/o~f At the LHC ener- 
gies a change of the regime, connected with the transition from shadowing to 
antishadowing is expected. 
The properties of multiplicity distribution of secondaries P(n) at high values of n remain 
among the topical problems in high-energy physics. As pointed out in a series of recent 
papers [0J , the underlying dynamics behind these rare processes may be quite different from 
the bulk of events. 

Our knowledge about high-energy multiplicity distributions comes from the data col- 
lected at the ISR, SppS collider (UAl, UA2 and UA5 experiments) and the Tevatron 
collider (CDF and E735 experiments). In should be noticed that the recent results from the 
E335 Collaboration taken at the Tevatron [0] do not completely agree with those obtained 
by the UA5 Collaboration at comparable energies at the SppS collider || (see Fig. 1). 

Notice that the delicate features of *&(z) at very large multiplicities, near the large- z 
edge can be better seen if the variable z is used instead of n. 



* E-mail: jenk@gluk.org 
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FIG. 1. E735 results on charged particle multiplicity distributions in full phase space compared 
with UA5 results. Data from the two experiments which were taken at nearly the same energy are 
rescaled by the same factor. 

On the theoretical side, it became common [f§], [0, 0, to approximate the observed 
distributions by the convolution of two binomial distributions, accounting for the general 
"bell-like" shape of P(n) with the observed structures ("knee" and possible oscillations) 
superimposed. 

One of the hottest issues in this field is the dynamics of very high multiplicities (VHM) 
||, close to the kinematical limit imposed by the phase space. The VHM events are very rare, 
making up only about 10~ 7 of the total cross-sections at the LHC energy, which makes their 
experimental identification very difficult. An intriguing question is the possible existence of 
a cut-of in the VHM region, beyond z = n/ < n >^ 5, where < n > is the mean multiplicity. 
In our opinion, a better understanding of the underlying physics can be inferred only in a 



model involving both elastic and inelastic scattering, related by unitarity. Such an approach 
has been advocated in a series of papers ||, flOfl , [ |TTf , summarized in ref. |12| . 



After a brief summary of the main ideas behind this approach, we analyze the relation of 
the distribution of secondaries and the behavior of the elastic and total cross sections with 
the possible transition from shadowing to antishadowing [[T3[ . 

We show that the existence of a cut-off at high multiplicities in the distribution ^(z) is 
related to the validity of GS and KNO scaling. 

The basic idea of the geometrical approach to multiple production, used in the present 
paper, is that the number of the secondaries at a given impact parameter p, n(p, s) is 
proportional to the amount of the hadronic matter in the collision or the the overlap function 
G(p,s) 

<n(p.s)>=N(s)G(p,s), (1) 

where N(s) is related to mean multiplicity, not specified in this approach, and G(p, s) is the 
overlap function, related by unitarity to elastic scattering 

Imh(p,s) = \h(p,s)\ 2 + G(p,s). (2) 

where h(p, s) is the elastic amplitude in the impact parameter representation. Unitarity, a 
key issue in this approach, enters both in the definition of the elastic amplitude and of the 
inelastic one (the overlap function). 



In the u— matrix unitarization (see [0 and references therein) 



Gin = Z OC x. : i \Ki (3) 

1 + 2\su + \u\ 2 

where u is the elastic amplitude (input, or the "Born term"). 

We use a dipole (DP) model for the elastic scattering amplitude, exibiting geometrical 
features and fitting the data. After u-matrix unitarization, the elastic amplitude reads (see 

ED 

(4) 

2 

where u(y, s) = ige~ y , y = -£-rr, and L = Ins. 
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Remarkably, the ratio of the elastic to total cross sections in this model, 

(l + g)Hl + g) [) 

fixes the (energy-dependent) values of the parameter g. Typical values of g for several 
representative energies are quoted in |T^| . 

Rescattering corrections to Gi n (p,s) here will be accounted for phenomenologically ac- 



cording to the following prescription (see []12] and earlier reference therein). 



G in (p,s) = \S(p,s)\G in (p,s), (6) 

where S(p, s)) is the S is the elastic scattering matrix, related to the u matrix by 

Q( \ 1 + iu(p, s) 

S{p,s) = - — r. 7 

1 — iu[p, S) 

This procedure is not unique. For example, it allows the following generalization (see 



12 1 and earlier reference therein) 

G m (p,s) = \S( Pl sTG a (p,s), (8) 

where a is a parameter, varying between and 1. 
We assume 

<n(p,s)>=N(s)Gf n (p,s). (9) 
The moments are defined by (see JE3] and earlier references therein) 



r k, ^_ N k {s)5G m {p,s){GT a {p,s)) k <?p 
Mj SG in {p,s)d* P UUj 

Now we insert the expression for the DP with the u-matrix unitarization (4) into (10) to get 

kt \ ^ _ N k (s)(l + g) [9 dx (fl + x^a x 



< n is) >- 



g Jo (1 + x) 2 \ v 1 — x ! (1 + x 
The mean multiplicity < n(s) > is defined as 

N(s)(l + g) [9 xdx ( l+x\a N(s) 



(^) a X . (11) 
Vl-J (1 + x) 2 1 J 



For the distributions we have 

P(n) = -/ 7 r?6[n-N( ) 7 ^ . (13) 

V ; 5 Jo (1 + x) 2 V U-x^ (1+x) 2 ; V ; 

Integration in (13) gives 

^(2) =< n > Pin) = -— — } {- r , 

^ K ' V ' g z{\ +x)[(l -x) 2 + 2axY 

where z — n/ < n > . 

Since the above integral is non-zero only when the argument of the 5 function vanishes, 

1 + x \ x 



n = N 

one gets a remarkable relation 



1 — x } (1 + a;) 2 ' 



(14) 



(l + x) 2 - a (l-x) a ' 

To calculate the distribution ^f?(z) one needs the solution of equation (14). It can be 
found explicitly for two extreme cases, namely a = and a — 1. Otherwise, it can be 
calculated numerically. 

The maximal value of z, corresponding to x = g (x varies between and g), can be found 

as: 

_ ag 

Zmax (i + s) 2 -1(i-s)K 1 j 

It can be seen from (15) that z max is a constant if g is energy independent. The ex- 
perimentally observed ratio <J e i/a t varies between 53 Gev and 900 Gev from 0.174 to 0.225, 
implying the variation of g from 0.489 to 0.702, uniquely determined by the above ratio. 
This monotonic increase of g(s) in its turn pushes z max (s) outwards, terminating when g 
reaches unity (according to [I2[ this will happen around 10 TeV, i.e. at the future LHC), 
whenafter the term |1 — g\ in (15) will start rising again, pulling z max (s) back to smaller 
values. I.e. z max (s) has its own maximum in s at g — 1. 

The unusual behavior of z max (s) is not the only interesting feature of the present ap- 
proach. This effect can be related to the behavior of the ratio a e i/cr t . As argued by Troshin 



and Tyurin (see |U|), cr e i/a t may pass the so-called black disc limit and continue rising in a 
new, " antishadowing" mode of the w-matrix unitarity approach (multiplicity distributions 
were not considered in that paper). According to the recent calculations |14| the transition 
from shadowing to antishadowing will also occur in the LHC energy region. 

To summarize, we found a regularity connecting the geometrical properties in high-energy 
dynamics (GS and KNO scaling) with the dynamics of the high-multiplicity processes. We 
showed, in particular, that exact geometrical, or KNO scaling, implying constant g in our 
model, results in a cut-off at large z of the distribution function ^l{z). Any departure from 
scaling (energy dependence of g in our model) shifts the point z m according to eq. (15). 
Within the present accelerator energy domain (ISR, SPS, Tevatron) g varies from about 0.5 
to about 0.8. It will reach the critical value g = 1 at LHC, where we predict a change of 
the regime: z max (s) will start decreasing and the black disc limit will be passed (which, as 
shown in and |14}| , is not equivalent to the violation of the unitarity limit, but means 



passage from shadowing to antishadowing [TJj and fl4]j). 

Finally, it should be noted that we use many model assumptions, decreasing the pre- 
dictive power of our calculations. These assumptions concern mostly the way absorption 
corrections are introduced and the assumption of the local (5 function) dependence of mul- 
tiplicities on the impact parameter. Both assumptions, as well as others can be modified. 
As a result we quantitive rather than qualitative changes in the results. 
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